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1 Introduction 

It is our common household experience that when the voltage drop across a 
fuse exceeds a limit, the fuse burns out. A fuse is nothing but a conductor 
which conducts uniform current under an applied voltage upto a certain limit 
beyond which it burns out and becomes nonconducting. This is called fuse 
failure. Similarly in a dielectric breakdown, a dielectric starts to conduct elec- 
tricity when the voltage drop across it attains certain threshold value. The 
above two phenomena are examples of breakdown process which is described 
broadly as the failure of a physical attribute when the perturbing force driving 
it goes beyond a limiting value. The most common example of the process is 
the breaking of a material at a high stress beyond its strength. 

Naturally occuring solids are almost always inhomogeneous and have de- 
fects like vacancies, microcracks or impurities. These defects are weak points 
across which stress fields or electric fields or current densities concentrate. 
Depending on the geometry of the defect, the concentrated field can be very 
high. If the field or density exceeds the material fatigue limit, failure nucle- 
ates locally around the defect and starts to propagate. The propagation can 
be arrested, or can spread to the entire system leading to the global failure 
of the system, depending on the field strength and the defect structure in 
the material. The failure of a material, thus, depends largely on the disorder 
present in it. 

The simplest form of the weak points is substitutional disorder which can 
be realised by the inclusion of nonconducting material in a conducting material 
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or vice versa. Since we are interested in the macroscopic properties of global 
failure, we consider failures at length scale larger than the regions in which 
the local failure appears (i.e., atomic distances). At this scale, the disordered 
fuse system can be modelled by a lattice whose bonds are conducting with a 
certain probability p and nonconducting otherwise. 

For such models, one can apply the principles of percolation theory [45]. 
Below the percolation threshold p c (for p < p c ) the system is not connected 
globally by channels of conducting bonds and so nonconducting, whereas for 
p > p c the system is conducting (at least one continuous path exists across 
the system via the conducting bonds). 
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Fig. 1. Phase diagram of a mixture containing p fraction of conductors and (1 — p) 
fraction of insulators at random. On the left side, the sample is insulating for voltage 
V less than the breakdown voltage 14 (p) and conducting otherwise. On the right 
side, the system is conducting for current I less than the fuse current If(p) and 
insulating otherwise. 

For p > p c , when conducting channels span across the system and the 
system is conducting, one can ask what is the fuse current If required so that 
the system becomes nonconducting. At I > If{p) not a single spanning path 
of conducting bonds exists. At p = 1 all the conducting bonds are present 
and //(l) normalized by the sample size is simply the fuse threshold current 
for each of the conducting bonds. For p < p c , the conducting bonds do not 
form a continuous path across the sample, and the sample is insulating. With 
increasing voltage V across the sample, one can get a continuous path through 
original conductors and broken dielectric for V > Vb(p), where Vb(p) is the 
dielectric breakdown voltage. At p = 0, all the non conducting bonds are 
there in the system and Vf,(0) normalized by the system size is simply the 
breakdown threshold of each of the nonconducting bonds. On the other hand, 
as p — > Pc, both the dielectric breakdown voltage and fuse failure current tend 
to zero. The corresponding phase diagram is shown in fig. 1. 

In this article, we will review the basic ideas in breakdown problems in 
terms of the fuse failure and dielectric breakdown problems and then discuss 
how these ideas can be extended to the breakdown in quantum mechanical 
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systems such as Anderson insulators. In disordered electronic systems above 
two dimensions, the electronic states below the mobility edge are all localised 
(the system with Fermi level within this range behaves as an insulator) and the 
states above the mobility edge are extended (the system turns to a conductor 
for Fermi level in this range). The Anderson transition from insulating to 
conducting phase across the mobility edge is well studied [26]. We discuss 
here the possibility of a breakdown from insulating to conducting phase by 
applying strong electric field and compare the quantum breakdown with the 
classical dielectric or fuse breakdown in disordered materials. 



2 Analysis of the fuse problem 
2.1 Disordered fuse network 

In a pure conductor placed between two electrodes with a potential difference, 
the field lines within the conductor are all parallal to each other and perpen- 
dicular to the electrode surfaces. In presence of a disorder in the form of an 
insulating region, the field lines get deformed around the defect. As a result 
in the vicinity of the defect current density increases to i e from i, the current 
density value far away from the defect. So one can write 



i e =i{l + k) (1) 

where k is the enhancement factor which depends on the geometry of the 
defect. As an example, k becomes l/b for an elliptic defect of semi-major axis 
I and semi-minor axis b [12]. For electrode surface area S total current I is 

Failure occurs for the first time when i e becomes equal to i$ which is the 
fatigue limit of the sample material. The failure current is then 

'/ = T~ TT" <3) 



Larger the enhancement factor k, smaller is the failure current If. For an 
example, if I >> b for an elliptic defect, If may get reduced by a large extent. 
So the presence of defects in the material facilitates failure and the presence 
of sharper edges of defects make the system more vulnerable to failure. The 
failure makes the defect bigger and hence k larger. This means that the current 
density around the defect enhances further causing failure again. The process 
causes rapid failure of the whole sample. This means that external voltage for 
the global failure is the same as the voltage for the first local failure and a 
local failures once started leads to the failure of the entire sample. This type 
of failures are called brittle failures . 
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So far we have discussed the influence of a single defect of a regular size 
within the sample. Natural and engineering samples usually contain a large 
number of defects of irregular shapes and sizes at random. To get some quanti- 
tative estimate of the failure criteria in terms of defect parameters we idealize 
the solid to lattice model and use percolation theory for defects [37, 36, 18]. 

We start with a hypercubic lattice in 2d or 3d with all conducting bonds. 
The simplest defect can be introduced by removing one bond parallal to the 
direction of the current flow. In this case, failure current is calculated as 



in 2d, where Li is the fracture current of the lattice in absence of any defects 
and L is the linear size of the lattice. The enhancement factor is 4/ir here. We 
can introduce randomness in shape and size of the defect by removing (1 — p) 
fraction of the bonds randomly. It is no longer possible to determine the most 
vulnerable defect and to calculate the enhancement factor. We consider two 
limits: (1) dilute limit (p — > 1) when the defect density is small and (2) near 
the critical point p c where the defect density is large and beyond which the 
lattice loses its connectivity. 

Dilute limit (p — > 1) 

In this limit there are only a few isolated defects (insulating bonds) in the 
sample. The current density around a defect is effectively independent of the 
presence of other defects. Under such considerations the most vulnerable of 
the defects will lead to the failure of the system. So our primary task is to 
identify the most probable dangerous defect (the weakest point which causes 
the largest concentration of the current density). 

The ensemble of n successively removed bonds, far from boundaries [18, 
19], in a plane perpendicular to the current flow, act as dangerous defect. In 
2d it is a linear defect and in 3d the defect has the shape of a shape of disc. 
The current through the conducting bonds at the immediate neighbourhood 
of the defect is 

i e = i(l + k^n) (in 2d), i e = i(l + k 3 n^) (in 3d). (5) 

Here i is current density through the bonds far away from the defect. En- 
hancement factor k contains y/n term because current is diverted by n defects 
to spread uniformly around the perimeter of the disc, which is proportional 
to \fn. The probability of appearance of a defect of n successive insulating 
bonds is 



where L is the lattice size and d is the dimension of the lattice. L d provides 
the number of places that the defects can occupy. The approach of P(n) w 1 
gives the size of most probable dangerous defect: 



(4) 



P(n)~(l-p) n L d , 



(6) 
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■InL (in 2d), 



In(l-p) 

Combining (5) and (7) i e becomes 

l + fc 2 

= t, 1 + fc 3 



In(l-p) 



InL (in 3d). (7) 



2 InL 
ln(l-p) 

3 InL 

"in(l-p) 



(in 2d), 
(in 3d). 



(8) 



Here, the total current is iL^ d x \ Equating i e with the threshold value io, the 
expression for the failure current becomes 



If = 



i L 



l + 2fc 2 [ 



In L 



(in 2d), 



Lln(l-p)J 

J ^ = i i 7^ r In L ,■ ( m3rf )' 



(9) 



From the equations it is clear that as p — > 1,7/ reduces to i^L^ 1 , the value 
of the current in the lattice with all the bonds conducting. The slope of the 
If(p) vs. p curve at p = 1 is infinite (see fig. 1). It is expected because even 
the presence of a single defect enhances the process of breakdown through 
cascade. The most important thing is the size dependence of //. For large 
enough L and p not very close to 1 (such that the absolute value of ln(l — p) is 
not too large) the failure current per bond if(= I//L d_1 ) decreases as 1/lnL 
and l/^lnL) 1 / 2 in 2d and 3d respectively. 



Critical behavior (p — > p c ) 

Near criticality (p — > p c ), the material is strongly disordered and can be de- 
scribed by node-link-blob picture [43, 41, 17] of percolation theory (see fig. 2). 
Close to p c the conducting part of the material extends over the sample to 
form percolating cluster which is self-similar up to the length scale £ p , the 
percolation correlation length. The geometry of the cluster at length scale 
£ p appears same as that of the original cluster at smaller scale. This means 
that the infinite cluster may be divided into cells of size £ p . Each cell consists 
of backbone bonds (which takes part in current conduction) and numerous 
dangling bonds (bonds which do not take part in current conduction). Only 
backbone is important because it takes part in conductivity. The backbone is 
made of two kinds of bonds: multiply connected bonds appear as blobs; and 
singly connected bonds (called red bonds), appear as links. Distribution of 
current in the sample is solely determined by links; being singly connected, 
each link has to carry the full current inside a cell. Since at threshold (p c ) 
correlation length (£ p ) spans the sample (large enough), it is reasonable to 
assume that the failure current approaches zero as p — > p c . 




Fig. 2. A portion of the node-link-blob superlattice model near p c . The distance 
between two nodes of the lattice is £ p , while chemical length of the tortuous link of 
the super lattice is L c . 



To determine quantitatively the critical behavior of the failure current If , 
we consider the voltage V applied across the node-link-blob network. Since 
V is distributed among L/£ p number of links (or cells) in series, the average 
voltage across each link is Vl ~ £, P V. The resistance of the sample R is 
related to average link resistance Rl through R ~ ^~ 2 Rl, since there are 
L^ 1 number of links in series in the length L of the sample and there are 
^p 1 number of parallal such links. So the mean current in a link appears to 
be il ~ £ p l ~ 1 V/R through the relation il = Vl/Rl- V equals to failure or 
fuse voltage Vf when reaches the maximum current which a link can stand. 
One can assume fuse voltage behaves as 

V f ~{p-p c )-*'. (10) 

Using the relation £ p <~ (p — p c )~ u and R <~ (p — p c )~ tc , tf comes out as 
tf = t c — (d — l)v. From the relation If = Vf/R, we get 

If-ip-p^-V". (11) 

The values of the correlation length exponent v are 1.33 in d — 2 and 0.88 
in d = 3 and the values of the conductivity exponent t c are 1.33 in d — 2 and 
2 in d = 3 [45] . So it is clear that fuse voltage Vf attains finite value in 2d and 
diverges with the exponent 0.2 in 3d, in contrary to the failure current which 
always approaches zero as p — > p c . 

Influence of the sample size 

Size dependence of // is related to the notion of the most dangerous defect 
present in the sample. In the present context the most dangerous defect is 
a 'cell' of the infinite cluster with length £ p in the direction parallal to the 
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applied voltage and l max in the perpendicular direction. The total probability 
of having a defect of size Z is, P — g(l)(L / 'S, p ) d , where g(l) is the probability 
density of defect cluster of linear size I. Percolation theory predicts [45, 38] 
that 

p ~ -'H) (i) d - (12) 

lmax is obtained when P«l and 

lmax ~ £,p hl£. (13) 

Now the current that flows through the side link of the defect is proportional 
to (I max) d 1 I and one obtains 

I f „ [P Pc > (14) 

So this is the correction over equation (11) due to finite size of the sample. 
Bergman and Stroud [7] gave an idea about the competition between extreme 
statistics and percolation statistics. The extreme statistics (size dependence 
of the most probable defect and failure current) is expected to dominate for 
In L > £ p (or, when Lifshitz scale is greater than the connectedness correlation 
length). The dominance of extreme statistics is expected for p far away from 
percolation threshold when the correlation length £ p is small. 

Li and Duxbury [27] proposed the dependence of If on L through (In L)~^ f . 
The approximate range of ipf is 

<V/<1- (15) 



2(d- 1) 

Considering the results of the dilute limit and in the critical region, the com- 
bined form of (9) and (11) is 



If = In 





(P-Pc) 


4> f 




_(1-Pc)_ 





1 + K 



In(l-p) 



(16) 



The value of the different exponents (see table 1) and the constant K depend 
on the dimension and on the type of percolation. From (9) and (14) one can 
see that the combined result is valid only for 2d. The expression has three 
obvious features: 

1. For p = 1, If = Icj, as expected. 

2. Near p = 1, (p — p c ) is almost constant and we get back the expression (9). 

3. Near p c , the denominator of (16) is of the order of unity and we recover 
the expresion (11) with <j>f = (d — l)u. 
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2.2 Distribution of the failure current 
dilute limit (p — > 1) 

In random fuse networks the failure current If shows large sample to sample 
fluctuations. Since the failure current is determined by the weakest defect in 
the sample, the fluctuations in the failure currents do not come down with 
the system size. The average failure current is not a self-averaging quantity. 

Distribution of the failure currents in a system of size L follows as (see ref. 
[12]) 



Pi, (I) = 1-exp 



-A d P d cxp < -cL41nP 



Io 



1 



d-1 ' 



(17) 



The derivative of this cumulative failure probability distribution F L (I) 
with respect to current / provides the current / at which the system fails. 
Certainly, at most probable failure current the derivative becomes maximum. 
Here io is the failure current for pure sample. It may become obvious with a 
simple calculation that If appears to be the most probable failure current (as 
was assumed in calculation) only in the limit of large enough system size. 

Though the current / can vary from zero to infinity, the form of Fl (I) is 
meaningful only for / upto Iq- One should expect the value unity for Fl{I) 
at I = Io and I = oo, but it is true only for large system size. Fl(I) suffers 
from size and defect concentration dependence through If (see equation 9). 

The above expression is reffered to as the Gumbel distribution [23]. An- 
other well known one which is very often used in engineering is Weibull dis- 
tribution 

I 



F L {I) = 1-exp 



-rL a 



(18) 



Here m is a constant and for large m (say more than 5) If refers to most 
probable failure current. 

At critical region (p — > p c ) 

Near criticality the cumulative failure distribution function is 



F L (I) = 1-exp 



-A' d L d exp - 



k' (p- Pc y 

1 

J d-1 



(19) 



where A' d and k' are two constants. 



2.3 Continuum model 



One can extend the ideas of lattice percolation to continuum conducting 
medium. A material at the scale of the size of the defects can be looked 
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upon as a continuous field with some defects as inclusions in the field. In con- 
tinuum model of lattice percolation, insulating spherical holes (circular holes 
in 2d) are punched at random as defects in a uniform conducting sample. The 
holes can overlap (Swiss-cheese model) and two non-overlapping neighbour- 
ing holes have a conducting region between them. These regions constitute 
conducting channels of cross-section S. Somewhat similar to the breakdown 
problem discussed above, the transport properties of any such channel de- 
pend on the transport capacities (cross-section and length) of the narrowest 
part (the weakest bond) of the channel. With some reasonable assumptions, 
one [24] can express the transport capacities of the weakest region in terms 
of percolation cluster statistics on the lattice (particularly in terms of the 
percolation correlation length £ p ). 

Both the discrete and continuum type are almost same in the dilute limit 
(p — > 1). So all the results derived earlier for discrete model are valid here. But 
near criticality an infinite percolation cluster, with the links of mean length 
£ p and of different cross-sectional width 5, is formed. The backbone [17, 41] 
of this cluster is represented by a superlattice (see sec. 2.1) of tortuous singly 
connected links and blobs crossing at nodes at a separation 

Zv~\p-Pc\-". (20) 

The chemical length between any two nodes is 

L c ~\p- Pc \-t. (21) 

For singly connected bonds (or sites) on the percolating backbone, Q = 1 
in all dimensions [24] except d = 1 [45, 15]. There arc (L/^ p ) d ~ 1 number of 
parallal links exist between two electrodes, where L is the size of the sample. 
The current in a link is given by %l = [£,p/ L) d The current density i in a 
channel of cross-section 5 is given by 

cd-lj 

* - frr (22) 

The maximum current density is obtained for minimum cross-section 8 m in 
which is inversely proportional to the shortest chemical length L c [24, 4, 6]. 
If ig is the threshold current density at which the material fails, then 

io-^- 1 ^- 1 // (23) 

and we get, 



If~(p-Pc) v+1 (in 2d), J/~(p-Pc) 2( " +1) (in 3d). (24) 

Thus the exponents for the failure current are higher than those of the 
discrete model. 
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Using If — Vf/R and R ~ \p — p c \~ tc with equations (20) and (21) the 
expression for fuse voltage becomes Vf ~ \p— p c \ tf with tf equal to (d— 1)(^ + 
1) — t c [12]. i c is the conductivity exponent: t c ~ 1.3 and 2.5 in 2d and 3d 
respectively. Consequently, £/ becomes approximately 1 and 1.3 in 2d and 3d 
respectively. Unlike discrete percolation where failure current at p c attains a 
finite value in 2d and vanishes only in 3d, in continuum percolation failure 
voltage always vanishes at p c . 



Table 1. Theoretical estimates for the fuse failure exponent <f>f 



Dimension 


Lattice percolation 


Continuum percolation 


2 


„(= 4/3) 


i/ + l(=7/3) 


3 


2i/(~ 1.76) 


2(i/ + l)(~ 3.76) 



2.4 Electromigration 

Electromigration is an example where we find the practical application of 
the concepts we have hitherto discussed. Miniaturization of the circuits and 
gadgets have become the norm of the day. Very thin metallic films are used 
as interconnects among the active parts of devices, resulting in large current 
densities through the metallic films. Having large enough momentum, the 
free electrons become able to displace metallic ions from their equillibrium 
positions. Thus depending upon the material a net material transport occours 
[48] through grain boundary diffusion, surface diffusion or lattice diffusion. 
This ionic displacement and the accumulated effect of material transport due 
to large current densities (more than 10 4 A/cm 2 ) are called electromigration. 
Over time such electromigration leads to void formation at the cathode and 
extrusion at the anode in thin film interconnects. Such void formation and 
hillock formation cause an open circuit and a short circuit to neighbouring 
connecting wires respectively. The problem of material accumulation can be 
suppressed by the layers of other material around and above the interconnects. 
The problem of open circuit due to void formation has received much attention 
as the random failures of several interconnects finally lead to failure of the 
entire device. 

Let a fraction (1 — p) of the resistors are removed from random resistor 
network (say in 2d). A random walker starts its walk from one of the lateral 
sides (those without electrodes) of the network, and jumps from one cell to a 
nearest neighbour cell by crossing the bonds irrespective of occupied or unoc- 
cupied bonds, intending to reach the opposite side, with the constraint that it 
never visits a cell more than once. Let the whole path consists of no occupied 
and n\ unoccupied bonds. For a given configuration of missing resistors the 
minimum value of no is the shortest path.The mean shortest path < no > (for 
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theoretical study see ref. [14, 46]) can be realised by considering the average 
over a large number of configurations. An electromigration induced failure of 
the network can be realised if the walker fuses the resistors as it crosses them. 
The criterion for a resistor to fuse is not a particular value of current (as in 
the fuse model) , rather a particular value of thershold charge Qo that attains 
the resistor from the time of application of the constant current Io to the 
network (for details of the model see ref. [10, 52]). If t\ is the failure time of 
an arbitrary resistor of the network which is subjected to a constant current 
Jo, then the resistor must satisfy 

/ 1 1(t)dt = Qo (25) 
Jo 

to fail. 

So shortest path for a given impurity concentration is the path which 
corresponds to the smallest number of resistors to fuse. The problem is to 
study the variation of failure time of the whole network r with p. For a network 
having an isolated defect of length n, r is given by 

T = (L-n)Q /I , (26) 

where LQq/Iq is the r value for the pure network. In pure limit the mean 
failure time < r > is related to the longest probable defect with size n c . This 
means that the bonds which fuse are also the bonds with largest current. 

r decreases as p decreases from pure limit, and approaches zero at p c , since 
the number of bonds to be broken goes also to zero. The shortest path no, 
the failure current If, the number of broken bonds Nf and the failure time 
< r >, all these quantities tend to go to zero at p c with the correlation length 
exponent v. 

2.5 Numerical simulations of random fuse network 

The numerical simulation of failure of random fuse network in 2d was carried 
out by de Arcangelis et al. [16]. Here, one starts with a lattice, the bonds of 
which are conductors and present with a probability p. An increasing external 
voltage is applied across the lattice and the voltage Vj when there is first local 
failure is recorded. The fused bond is removed and the voltage is raised till 
there is a second failure at Vj. The process is continued till the global failure 

at a voltage vj m occurs. The variation of Vj and vj m with p are studied. 
Vj decreases with decreasing p till p ~ 0.7 whereupon it attains minimum 
value and then it starts increasing again. Vj m , on the contrary, increases 
monotonically and behaves almost identically as Vj for \p — p c \ < 0.08. Both 
approaches p c with the diverging exponent 0.48. It seems that Vj exhibits 
pseudodivergence in 2d. 
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Duxbury et al. [19] performed similar type of simulations covering the 
whole range of p from 1 to p c — 0.5. Their simulation results concerning If 
vs. p graph fit well with the iterpolation formula (16) with the exponent value 
4>f = 1, whereas the theoretically predicted value is 1.33. This may be due to 
the finite size effect. 



3 




10 



30 
L 



Fig. 3. L/If vs. In!/ expressing their linear dependence (after Duxbury et al. 
[19]). The curves from top to bottom correspond to initial impurity probability 
p — 0.6, 0.7, 0.8 and 0.9 respectively. 



They also demonstrated the finite value of Vj at p c by looking into the 
variation of If and the conductance near p c . Following equation (9) they 
checked the linear dependence of L/If on InL successfully for several values 
of p varying L from 10 to 200 and determined ijif = 1 from their slopes. The 
slope increases as p approaches p c . Instead of (17) they preferred 



F L (V) = l-exp 



-AL 2 exp 



KL 



(27) 



as cumulative distribution function of failure voltage. 

Arcangelis et al. [16] and Duxbury et al. [19] determined the number of 
fused bonds Nf upto complete failure and they found that Nf goes to zero 
algebraically as p — > p c with an exponent seemingly equal to the correlation 
length exponent v. 



3 Dielectric breakdown problem 



In dielectric breakdown problem the sample is insulating and conducting ma- 
terial acts as defects in the composite. The volume fraction (p) of the conduct- 
ing material is less than the critical volume fraction (p c ) so that the system 
overall is not conducting. The dielectric portions can withstand electric field 
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upto e c , at and beyond which they become conducting (local breakdown). 
Global breakdown occurs when the conducting portions span the sample with 
the succession of local breakdown under the influence of an increasing external 
field. 

In 2c? the solution for dielectric problem can be obtained from the solution 
of fuse problem with the aid of duality relation [33, 9]. We follow Bowman and 
Stroud [9] and consider the case where p is smaller than p c and the sample is 
macroscopically insulating. The equations for the induction vector D and the 
field E are: 



V-D = 

V x E = 
D(r) = e(r)E(r), 



(28) 



where e(r) is local dielectric constant of the insulating portion and E is irro- 
tational . With E = — v <j> equations (28) yields 



d_ 

dx 



e(r) 



dx 



+ 



d_ 

dy 



e(r) 



dy 



= 0, 



(29) 



where <p is the scalar potential. Now consider the dual composite of the origi- 
nal, where insulator phase of the original is replaced by conducting phase and 
vice versa. The sample is macroscopically conducting and the relevant quan- 
tities are current density i and electric field E. These satisfy the equations 



yi = o 

V x E = 
i(r) = <7(r)E(r) 



(30) 



Since i is divergenceless, i can be expressed as a curl of a vector potential V. 
V (V = V z = ip(x, y)) is chosen in such a way, only z-component of i vanishes. 
Assuming local conductivity as a(r) — l/e(r) equations (30) yields 



d_ 

dx 



C (r) 



dtp 
dx 



d_ 

dy 



C (r) 



chp_ 
dy 



0. 



(31) 



Comparing (29) and (31), one has 



(32) 



dip d<j> dtp def) 
dx dx ' dy dy 
With all these, the components of current density i in dual composite become 



dtp/dy = E y , 



-dtp/dx = —E x . 



(33) 



Thus we see that the magnitude of the current density in the dual compos- 
ite is equal to the magnitude of the electric field in the original one and the 
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direction is rotated by 90°. It can also be shown that the field E of the dual 
problem is equal to vector D of the dielectric problem and is rotated by 90°. 
Now the physical correspondence of the two pictures can easily be established 
[12,39]. 

All the results of fuse problem can be utilized for obtaining the solution 
for dielectric problem. For example, the equations (16) and (19) can be used 
in a straight forward way on replacing la and // by Vq and V& respectively 
and also (1 — p) by p. 



3.1 Dielectric breakdown problem 

As in the fuse problem, dielectric breakdown problem can be viewed on a 
discrete lattice with disorder in the framework of percolation theory. Basically, 
a lattice of insulating bonds is considered out of which p (p < p c ) fraction 
of bonds are conducting at random. For this purpose resistor (capacitor) of 
different resistivity (capacitance) [8, 9] ([5]) can be used. Resistors (capacitors) 
of smaller (higher) resistance (capacitance) are used for conducting bonds. 



Dilute limit (p 0) 

The problem is very similar to the fuse problem. The enhanced local field due 
to the presence of long defect (made of n number of larger valued capacitors 
(conducting bonds) perpendicular to the electrode, is (following Beale and 
Duxbury [5]) 

E e = E(l + kn), (34) 

where E is the externally applied field and k is the enhancement factor. It 
is valid for all dimension The probability to find a long defect made of n 
conducting bonds is 



P(n) 



p n L d . 



(35) 



The most probable defect size is 



lnp 



InL. 



(36) 



The enhanced field near the insulating bond adjacent to the most probable 
defect is now 

lnLV 



E P =E 



1+K d 



lnp J 



(37) 



When enhanced field E e of any bond attains threshold value of dielectric 
breakdown, local breakdown takes place and the bond becomes conducting. 
So breakdown voltage becomes 

E L 



1 



(38) 



where E is the breakdown field without defect. 
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Close to critical point (p — > p c ) 

In this limit, consider the conducting defects which are on the average 
the percolation correlation length £ p distance apart. Hence, average field is 
Vi/£ p = V/L, where V and V\ are externally applied voltage and potential 
difference between any two conducting defects respectively. The maximum 
attainable field is V\/a, where a is the minimum available separation (bond 
length) between the defects in the network. When the local electric field Vi/a 
reaches the bond-threshold value e c , breakdown occurs. Now the required av- 
erage electric field for local breakdown is E b — (a/£ p )e c . Near p c , £ p diverges 
as £ p (p) oc a(p — p c )~ v ■ So the critical behavior of the average breakdown field 
as derived by Lobb et al. [28] appears as 

E b ~{p c -p) v . (39) 

in all dimensions. This can also be derived [46] from the concept of minimum- 
gap g(p), as Eb can be considered as proportional (see sec. 2.4) to g(jp) and 
g(p) is proportional to t;^ 1 (see relation 20), the result above follows. Since 
Eb is inversely proportional to the linear dimension of the most vulnerable 
defect and the typical size of the defect or conducting clusters diverges at p c , 
Eb approaches zero in this limit. 



Influence of sample size 

We consider a critical defect which is a pair of longitudinal very closely spaced 
clusters of conducting material (almost linear in shape) of the order of size I 
situated one after one, parallal to the direction of the applied field E in a L d 
lattice. The field between these two clusters enhanced by a factor I time E. Far 
frompc these clusters appear with the probability (l/£ p )exp(—l/£ p ). Thus the 
most dangerous defect (the largest cluster) is of the order of l max ~ £ p \n(L d ). 
The first breakdown field El (field needed to break the first bond) scales as 
1 / lmax in a finite but large system [5] . Since, in this particular case breakdown 
field for the sample is same as for the first breakdown, Eb becomes 

E b ~ (40) 

This approximate expression for the average breakdown field was derived by 
Beale and Duxbury [5]. 

Larger sample has higher chance of having a larger vulnerable defect. So 
larger sample requires smaller electric field for breakdown and Eb becomes zero 
in the limit L — > oo. Breakdown path is always on the average perpendicular 
to the electrode. So the above result is independent of dimension. 

Due to the dominating behavior of (p c — p) v over InL, breakdown field 
tends to zero as p approaches p c . Very near to p c Chakrabarti and Bcnguigui 
([12] p. 66) proposed a scaling relation as 
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p (gc - vY , AU 

Eb ~ mW (41) 

According to Bergman and Stroud [7] , Eb may become size independent very 
close to p c , and a cross-over from extreme statistics to percolation dominated 
statistics may be there. 

Summary 

A general formula may be inferred for dielectric breakdown field as: 



Eb = E 



l + K 



— In p 



(42) 



with the exponent 4>b (=v for lattice percolation) dependent on the dimension 
and on the type of percolation. For the theoretically estimated value for 4>b 
see table 2. 



Table 2. Theoretical estimates for the dielectric breakdown exponent (f>b 



Dimension 


Lattice percolation 


Continuum percolation 


2 


u{= 4/3) 


i/ + l(=7/3) 


3 


0.88) 


+ 1.88) 



3.2 Distribution of breakdown field 

Distribution function Fl{E) in dielectric breakdown is defined as the proba- 
bility of breakdown of a dielectric system of size L in an external electric field 
E. Usually, Weibull distribution [50, 49] 

F L (E) = 1 - exp{-rL d E m ) (43) 

is used to fit the distribution function of breakdown and failure problems. Here 
r and m are constants. Duxbury and his co-workers [18, 19, 20, 5] argued that 
the distribution function for dielectric (and electrical) breakdown is given by 
Gumbell distribution: 



F L {E) = l-exp 



-AL d exp 



E 



(44) 



It can be derived [5] from a simple scaling argucment based on percolation 
cluster statistics [44, 21, 25]. It is valid for L » £ p or in dilute limit. The 
two expression (43) and (44) are qualitatively same if the Weibull exponent 
m is large. Though second one provides better fit to numerical simulation [5]. 
Sornette [42] argued that the expression (44) is not appropriate in continuum 
system having percolation type of disorder, rather, Weibull-like distribution 
is valid there. 
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3.3 Continuum model 

The model consists of dielectric material with the spherical (circular in 2d) 
conducting inclusions as defects with randomly positioned centres having the 
possibility to overlap. We assume breakdown field E b is proportional to the 
width 5. For breakdown numbers of link elements to be broken. With the 
same logic as in sec. 2.3, E b ~ t,,^ 1 ~ i" 1 ^" 1 or 

E b ~(p c - P r +1 (45) 

very near to p c [13, 28]. So continuum system is weaker than the discrete one. 
This is because the conductivities of the conducting channels increases as 
p — ► p C} whereas, conductivity of the bonds in discrete model is independent 
of p. 

3.4 Shortest path 

The concept of shortest path is already mentioned in section 2.4 in the context 
of electromigration. Just like before, we consider a walker starting its walk 
from one electrode with jump from one site to another and reaches the opposite 
electrode after executing a self-avoiding walk. The aim of the walker is to 
create a percolating conducting path by transforming a insulating bond to 
conducting bond when it jumps between two conducting sites separated by 
the insulating bond, no denotes the total number of insulating bonds which 
arc to be transformed into conducting bonds during walk across the sample. 
After completion of the walk, the sample experiences a continuous conducting 
path between the electrodes and in this way dielectric breakdown comes into 
existance. 

The shortest path for a given configuration is defined as the path with the 
smallest no-The normalised smallest path is given by g(p) =< no > /L, where 
< no > is obtained by considering the average of no over a large number of 
configurations. Some authors call g{p) by minimum gap. 

The behavior of g(p) has been studied [46, 20] extensively in 2d and 3d for 
regular and directed percolation scenario. As p increases from zero to p c , the 
minimum gap g{p) decreases from unity to zero with the correlation length 
exponent v. 

3.5 Numerical simulations in dielectric breakdown 

Stochastic models 

In these models, stochastic growth processes are considered which mimic 
the dielectric breakdown processes. For example, Sawada et al. [40] considered 
a random growth process where growth takes place in two ways: tips of the 
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pattern grow with a probability p and new tips (branching) appear with 
probability p n (here, p n < po). In their simulation the pattern appears as 
fractal and fractal dimension can be tuned by the parameter R (= p /p n ). 
However, results of such a simplification do not satisfy experimental results 
of dielectric breakdown [35] . 

To mimic the dielectric discharge pattern in gasses Niemeyer et al. [34] 
suggested a stochastic model where the breakdown pattern generated in turn 
determines the local electric field and the growth probability. The model could 
reproduce the fractal properties of dielectric breakdown process by numerical 
simulations. In this model the breakdown patttcrn starts growing from the 
centre of a lattice with insulating bonds. One electrode is placed at the centre 
and other one is placed at a long distance on the circumference of a circle. 
In one step only one interface bond (and a point) among all the nearest 
neighbors of the pattern breaks down depending upon the growth probability 
p and becomes the member of the pattern. The newly added bond becomes 
a conductor and the newly added site is shorted to the voltage of the central 
electrode. The growth starts from the centre and grows radially outwards. 
The growth probability p depends on the local field (potential) which in turn 
is controlled by the breakdown pattern via the relation 

p(i,3^i',3')= ^'K r (46) 

where the indices i,j and i/,jf represent the discrete lattice coordinates. The 
electric potential V is defined for all sites of the lattice by the discrete Laplace 
equation 

Vi,j = \{Vi+i,j + Vi- ld + V iJ+1 + Vij-t) (47) 

with the boundary condition V — for each point of the discharge pattern 
and V = 1 outside the external circle, rj describes the relation between local 
field and growth probability. The fractal structure of the pattern has been 
seen to obey 

N(r)~r d f, (48) 

where N(r) is the total number of discharge points inside a circle of radius r 
and df is the Hausdorff dimensions. In 2d one has df ~ 2, 1.89 ± 0.01, 1.75 ± 
0.02, 1.6 for i] = 0, 0.05, 1, 2 respectively. The structure tends to be more linear 
with larger i]. The observed value of d for i] = 1 is in good agreement with the 
experimental result (~ 1.7) (Niemeyer and Pinnekamp [35]) and the resulting 
figure is very similar to Lichtcnberg figure. 

The same pattern for 77 = 1 has been produced surprisingly by a differ- 
ent growth model: diffusion-limited aggregation model (DLA) of Witten and 
Sadner [51] (For review see Meakin [31]). Niemeyer et al. neither justifies the 
appearance of 77 nor provides any theoretical explanation regarding explicit 
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rule for breakdown. Many models (for details see ref. [39]) have been proposed 
but a model which fully describes dielectric breakdown in solids is still lacking. 

Deterministic models 

Contrary to Niemeyer et al. [34], Takayasu [47] introduced deterministic 
approach to produce dendritic fractal pattern (as is found in lightning) by con- 
sidering a priori spatial fluctuation on bond resistances (r^ = 8r; 8 G [0,1]) 
and the nonlinear irreversible characteristics of the resistance: once the po- 
tential difference across a resistance attain a pre-assigned threshold voltage 
v c , ri reduces to Sri (8 is a small positive quantity) and it's resistivity then 
never changes. This is the breakdown of a resistor in this model. Breakdown 
of a resistor induces successive breakdown of other resistors leading to the 
formation of percolation cluster of broken resistors. The pattern appears to 
be anisotropically fractal with the dimension df = 1.58 ± 0.12 in 2d. Due 
to anisotropy this df is less than the fractal dimension 1.89 of a percolating 
network in 2d. 

Family et al. [22] made the stochastic model of Niemeyer et al. [34] deter- 
ministic by attaching randomly a breakdown coefficient 9 (9 e [0, 1]) to each 
insulating bond. There are two versions of this model: In one model, at each 
time step an interface bond ij with the largest 8Vjj breaks down, whereas in 
the other model an interface bond breaks down with a probability 9V-j/p max , 
where r\ is an adjustable constant and p m ax is the largest value of 9V^ among 
all the interface bonds. The patterns appear as tunable (varying rj) fractal. In 
the first model stringy highly anisotropic (with branches) pattern with fractal 
dimension of about 1.2 appears in 2d, whereas in the second one pattern is 
strikingly lacking in anisotropy. The shape of the pattern is similar to those 
found in Niemeyer et al. [34] except the fractal dimension of 1.70 ± 0.05 (for 
r\ = 1) in 2d. The results for df appear same to that for DLA on a square 
lattice [30] and to that for the same model in a homogeneous medium (here, 
cross and square-shaped patterns appear) [22]. 

Direct simulations of dielectric breakdown from various groups confirmed 
the theoretical consideration for size and impurity dependence of breakdown 
voltages quite successfully. In these simulations, a discrete lattice of insulat- 
ing bonds is considered with p fraction of bonds as conductors. Each insulator 
breaks down to conductor at a threshold value of voltage drop v c . A macro- 
scopic voltage is applied and the voltage distribution throughout the lattice is 
computed by solving Laplace equation (47). The insulator with largest volt- 
age drop (at or above v c ) is converted to a conductor (incidence of first local 
breakdown). The voltage distribution is then recalculated, the second local 
breakdown is identified and the process continues. If at any step the applied 
voltage is not large enough to cause breakdown of any insulator, it is in- 
creased gradually. The simulation continues till a sample spanning cluster of 
conductors appears. The applied voltage Vj, at which global breakdown occurs 
divided by the sample size L is identified as breakdown field It has been 
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experienced that the breakdown field Ej, is same as the field E\ required for 
the occurrence of first local breakdown. 

Manna and Chakrabarti [29] determined the p dependence of both E^ and 
g(p) for the entire range of p below p c . They found (see figure 4) that both E\ 
and g(p) go to zero at p c with the exponent of value almost equal to 1. They 
argued that the exponent is actually v and the smallness of the value is due 
to the smallness of the lattice size (L = 25, where exact p c is not reachable). 
Bowman and Stroud [9] found that El vanishes at p c with an exponent equal 
to 1.1 ± 0.2 in 2d and 0.7 ± 0.2 in 3d for both site and bond problems. These 
are consistent with the correlation length exponent (v) values of 4/3 and 0.88 
in 2d and 3d respectivey. 




Fig. 4. The variation of first breakdown field Eb and the minimum gap g(p) with 
initial impurity probability p, (after Manna and Chakrabarti [29]). The two quan- 
tities behave same very near to the percolation threshold p c , whereas far from p c 
they behave differently. 



Beale and Duxbury [5] proposed a relation for first (local) breakdown field 

as 

Eh ~ A(p)+B(p) \nL (49) 
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Near p c , they found A(p) and B{p) vary like (p c — p)~ v (in the lattice of size 
L = 50, 70, 100) as is expected from eqn. (40). From the plot of ln{B(p) ln(p)} 
versus — ln(p c — p) their data shows v = 1.46±0.22 which is in good agreement 
with exact value 4/3 in 2d. Their data fits well to the Gumbel double expo- 
nential form (44) for cumulative failure distribution. Manna and Chakrabarti 
[29] and Beale and Duxbury [5] found the exponent of eqn. (42) to be about 
1.0 and 1.2 respectively. 

Acharyya and Chakrabarti [2] found that starting with a concentration 
p (p < p c ) of the conductors, the rate at which the insulating bonds breaks 
down to conductors as the electric field is raised in a dielectric breakdown 
diverges at breakdown voltage 14 . 
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Fig. 5. The plot of dielectric breakdown susceptibility vs. applied voltage for sev- 
eral sample sizes (after Acharyya and Chakrabarti [1]). The inset shows that the 
difference of the maximum susceptibility field and the minimum sample spanning 
cluster field decreases with increasing sample size 



This indicates that global breakdown process is highly correlated very 
near to V&. They defined a new quantity named as breakdown susceptibility 
X (= dn/dV). Here n(V) is the average number of broken bonds at an external 
voltage V (< Vb) for a fixed value of p. It is obvious that for a sufficient volt- 
age V , n(V) saturates to L d . \ exhibits a maxima at V h w hich is different 
from V h , , where V/ m is the voltage required to creat the last member of the 
sample spanning conducting cluster. V b approaches V/ m with increasing 
sample size L. It seems that x is a divergent quantity for infinite sample. So 
there exists a possibility to predict v/ m without going to the complete break- 
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down point of the sample. Statistics of the growing clusters upto VJ, (minimum 
external voltage required to have the global connection across the sample via 
the conducting sites) in 2d site percolation was studied by Acharyya et al. [3] . 
Vb is identified not only as the point of global breakdown but also the point 
of divergence of the rate of the various statistical quantities like total number 
of conducting sites, the average size of the conducting cluster and the number 
of such clusters. 



4 Zener breakdown in Anderson insulators 

From the microscopic point of view, current conduction requires the mobility 
of the electrons across the material and for composites of metals and dielectrics 
or insulators this means that there must be percolating paths of conducting 
material in the system (as long as electron is considered a classical particle). 
If p denotes the concentration of the conducting material in the composite 
and p c the percolation threshold, then from our previous discussion we know 
that the conductivity will go to zero at p c as \p — p c \ tc [45]. For p below p c , 
the composite is ceases to be conducting and the classical breakdown voltage 
needed to make the system conducting goes as \p — p c \^ b (see sec. 3.1). 

Studies on Anderson transition (see ref. [26] for example) show that elec- 
tron as a quantum particle cannot diffuse even through the geometrically 
percolating path due to the coherent back scattering (of the wave function) 
from the random geometry of the cluster in a dimension less than three. Since 
all the states on any percolating lattice gets localized (exponentially) electrons 
do not diffuse through the disordered (percolating) lattice. 

In a three dimensional disordered systems there exist a mobility edge (a 
sharp energy level in conduction band) e c below which all the states are lo- 
calized and above which all the states are extended. Insulting or conducting 
phase appears depending upon the position of Fermi level e / below or above 
the mobility edge e c . The (Anderson) transition from the insulating phase 
(ef < e c ) to metallic phase (e/ > e c ) (where the electrons are quantum me- 
chanically percolating that is their wave functions extend all over) across the 
mobility edge in a percolating disordered solid, has already been well studied 
[26]. In metallic phase the conductivity increases as |e/ — e c \ tq ; t q ^ t c . 

In the case of a classical nonpercolating (p < p c ) system, an application 
of an external electric field Eb forces the system to undergo breakdown fol- 
lowing the relation Eb ~ \p c — p\ v (see sec. 3.1). In analogy with the above 
classical dielectric breakdown problem one can think of a problem of quantum 
dielectric breakdown (possibility of the appearance of a conducting path) in 
Anderson insulating (e/ < e c or quantum mechanically non-percolating) phase 
with an application of strong electric field. This may be thought as a kind of 
complementary problem of the Zener breakdown in insulator. 

In an insulator, states are localised exponentially within the atomic dis- 
tance a and effective band width w is in general larger than a. There exists a 
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possibility of tunneling of the states across the band gap if the energy gained 
by an electron in travelling an atomic distance a in an electric field E is larger 
or equal to the band gap Aef, (or when eEa > Ae b ). Here e denotes the 
electronic charge and the band gap is the separation between conduction and 
valance band. In Zener breakdown, bands get effectively tilted in strong elec- 
tric field E in the direction of the field, reducing the band gap effectively from 
Aeb to Aeb r . If the effective width of the reduced band gap (w = Acb r /eE) 
becomes of the order of atomic distance a interband tunnelling takes place 
for w < a and the insulation breaks down. So the Zener breakdown field Eb 
scales linearly with the band gap: Eb = Ae&/ea. 

Similar kind of breakdown of insulation can occur in the case of Anderson 
insulators with an application of strong electric field E in more than two 
dimension. States are localized (exponentially) with in the localization length 
£ q , where £ g varies as £ g ~ |Ae m | _l/,! . Here Ae m is the mobility gap; Ae m 
< , — Ef. One can now easily think of a possibility of tunnelling across the 
mobility gap, if the energy gained by an electron in travelling a distance £ g 
in field E is of the order of the mobility gap Ae m (or when eE^ q > Ae m ). 
Thus in contrary to the standered Zener breakdown in semiconductors, here 
breakdown field E b scales almost quadratically with mobility gap as [11] 

E b = ~ |Ae m | T - T g = l + v q (50) 

{v q ~ 0.9 [26, 32], giving T q ~ 1.9 in 3d). 



ot stales 




Fig. 6. Schematic density of states for Anderson insulators (d > 2) shown in (a), 
where the Fermi label e/ is below the mobility edge e c (metallic phase for e/ > e c ). 
For strong electric field E, the band of (localized) states get tilted and tunnelling 
occurs when the effective width w(= \ Ae mr \/eE) of the mobility gap Ae mr , is less 
than or equal to the localization length £ q . 
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The tunnelling probability g(E) decreases exponentially with the width 
of the effective barrier as g(E) ~ exp(~w/^ q ) ~ exp(—Ac m /eEt; q ) ~ 
ea;p(— \Ae m \ Tq /E). The cumulative failure probability of Gumbcl form (same 
as for fracture [20, 18, 19] ) for a sample of size L under field E is given by 
[11] 



F L (E) ~ 1 - exp[-L d g(E)} 



1 — exp 



-L d exp 



E 



(51) 



This gives 



E b 



|Ae r 



InL 



(52) 



as the size dependence of the typical breakdown field for and above which 
Fl(E) is significant. 



5 Conclusions 

We have discussed the classical failure of the fuse systems, the dielectric break- 
down and the quantum breakdown in the Anderson insulators. We have dis- 
cussed how the extreme value statistics and the resulting Gumbel distribu- 
tion arises in breakdown and failure processes, especially when the disorder 
concentration is low. At high concentration of disorder near the percolation 
threshold, we have discussed how the cross-over might take place from extreme 
value to percolation statistics. We have discussed the system size dependence 
that arises at the distribution of the failure current at low disorder regime. 
Finally, the extension of Zener breakdown phenomenon for band insulators to 
the disorder-induced Anderson insulators has been discussed in sec. 4. 
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